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Abstract 

Inclusive KgKg production in ep collisions at HERA was studied with the 
ZEUS detector using an integrated luminosity of 0.5 fb" 1 . Enhancements in 
the mass spectrum were observed and are attributed to the production of 
/ 2 (1270)/ajj(1320), /a(1525) and / (1710). Masses and widths were obtained 
using a fit which takes into account theoretical predictions based on SU(3) sym- 
metry arguments, and are consistent with the PDG values. The /o(1710) state, 
which has a mass consistent with a glueball candidate, was observed with a sta- 
tistical significance of 5 standard deviations. However, if this state is the same 
as that seen in 77 — > KgKg, it is unlikely to be a pure glueball state. 
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1 Introduction 



The existence of glueballs is predicted by QCD. The lightest glueball is expected to have 
quantum numbers J PC = ++ and a mass in the range 1550-1750 MeV [1]. Thus, it can 
mix with qq states from the scalar meson nonet, which have 1 = and similar masses. 
Four states with J PC = 0++ and I = are established [1]: / (980), / (1370), / (1500) and 
/o(1710), but only two states can fit into the nonet. In the literature, the state /o(1710) is 
frequently considered to be a state with a possible glueball or tetraquark composition [2]. 
However, its partonic content has yet to be established. 

The ZEUS Collaboration previously observed [3] indications of two states, /2(1525) and 
/o(1710), decaying to K^Kg final states in inclusive deep inelastic scattering (DIS) events. 
The statistical significance of the observation did not exceed three standard deviations. 
The state in the 1700 MeV mass region had a mass consistent with that of the / (1710); 
however, its width was significantly narrower than that quoted by the Particle Data Group 
(PDG) [1]. 

The results presented here correspond to the full HERA luminosity of 0.5 fb" 1 and super- 
sede the earlier ZEUS results. The measurement of the Kg Kg final states is presented in 
a kinematic region of ep collisions dominated by photoproduction with exchanged photon 
virtuality, Q 2 , below 1 GeV 2 . The data allow the reconstruction of the K S K S final states 
with much larger statistics than previously used. 



2 Experimental set-up 

The data were collected between 1996 and 2007 at the electron-proton collider HERA 
using the ZEUS detector. During this period HERA operated with electrons or positron^] 
of energy E e = 27.5 GeV and protons initially with an energy of 820 GeV and, after 1997, 
with 920 GeV. 

A detailed description of the ZEUS detector can be found elsewhere [4] . Charged particles 
were tracked in the central tracking detector [5] which operated in a magnetic field of 
1.43 T provided by a thin superconducting solenoid. Before the 2004-2007 running period, 
the ZEUS tracking system was upgraded with a silicon Micro Vertex Detector (MVD) [6]. 
The high-resolution uranium-scintillator calorimeter (CAL) [7] consisted of three parts: 
the forward, the barrel and the rear calorimeters. 



Here and in the following, the term "electron" denotes generically both the electron (e ) and the 
positron (e + ). 
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3 Event sample 



A three-level trigger system [4, 8] was used to select events online. No explicit trigger 
requirement was applied for selecting K S K S events. The photoproduction sample is dom- 
inated by events triggered by a low jet transverse energy, E?, requirement (E? > 6 GeV). 
Deep inelastic scattering events were triggered by requiring an electron in the CAL. 

Events were selected offline by requiring | Z vtx |< 50cm, where Z vtx is the Z— coordinate 
of the primary vertex position determined from the tracks. The average energy of the 
total hadronic system, W, of the selected events was ~ 200 GeV. The data sample was 
dominated by photoproduction events with Q 2 < 1 GeV 2 . 

4 Reconstruction of Kg candidates 

The Kg mesons were identified via their charged-decay mode, K s — > 7r + 7r~. Both tracks 
from the same secondary decay vertex were assigned the mass of the charged pion and 
the invariant mass, M(tt + it~), of each track pair was calculated. The Kg candidates were 
selected by requiring: 

• M(e + e~) > 50MeV, where the electron mass was assigned to each track, to eliminate 
tracks from photon conversions; 

• M(p7r) > 1121 MeV, where the proton mass was assigned to the track with higher 
momentum, to eliminate A and A contamination to the Kg signal; 

• pr{Kg) > 0.25 GeV and ^(iCg)] < 1.6, where px{Kg) is the transverse momentum 
and rj(K s ) is the pseudorapidity; 

• O2D < 0.12 rad (6*30 < 0.24 rad), where 62D (#3L>) is the two (three) dimensional 
collinearity angle between the i^-candidate momentum vector and the vector defined 
by the interaction point and the Kg decay vertex. For 82D, the XY plane was used. 

The cuts on the collinearity angles significantly reduced the non-i^ background in the 
data during the 2004-2007 period. These cuts were necessary due to the extra material 
introduced by the MVD. After all these cuts, the decay length distribution of the resulting 
Kg candidates peaked at ~ 2 cm. 

Events with at least two Kg candidates were accepted for further analysis. More than two 
Kg were allowed in one event, unlike for the previously published result [3], and all distinct 
combinations of K S K S were used. In the mass range of 481 < M(tt + 7t~) < 515 MeV the 
number of K% candidates is 1258399. 
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Figure [T] shows the invariant-mass distribution of K s candidates. A fit over the whole 
mass range including a first-order polynomial was used to estimate the background con- 
tribution at ~ 8%. The central region was fitted with two bifurcated Gaussian functions 
to determine the mass and width of the Kg meson. For the HERA II data, corrections 
were applied to take into account the extra dead material introduced into the detector. 
After the corrections, the mass and width of the Kg were compatible with the PDG value 
and detector resolution, respectively. 



5 Results 

The KgKg invariant-mass distribution was reconstructed by combining two Kg candidates 
selected in the mass window 481 < M(n + 7r~) < 515 MeV. Tracks used for the K S K S pair 
reconstruction were required to be assigned uniquely to each Kg in the KgKg pair. 

Figure [2^ shows the measured K S K S invariant-mass spectrum. Three peaks are seen 
at around 1300, 1500 and 1700 MeV. No state heavier than the / (1710) was observed. 
The invariant-mass distribution, m, was fitted as a sum of resonance production and a 
smoothly varying background U(m). Each resonant amplitude, R, was given a relativistic 
Breit-Wigner form [9]: 

where Mr and Yr are the resonance mass and width, respectively. The background 
function used was 

U(m) = m A ■ exp (—Bm) , (2) 

where A and B are free parameters. The K S K S mass resolution is about 12 MeV for the 
mass region below 1800 MeV and its impact on the extracted widths is small compared 
to the expected widths of the states [3]. Therefore, resolution effects were ignored in the 
fit. 

Two types of fit, as performed for the reaction 77 — > K S K S by the L3 [10] and TASSO [9] 
Collaborations, respectively, were tried, using Eqs. (1) and (2). Fit 1 is an incoherent 
sum of three Breit-Wigner cross sections representing the /2(1270)/a2(1320), /2(1525) and 
/ (1710) plus background. Fit 2 is motivated by SU(3) predictions [11]. The decays of 
the tensor (J p = 2 + ) mesons /2(1270), a2(1320) and /2(1525) into the two pseudoscalar 
( J p = 0~) mesons K°K° are related by SU(3) symmetry with a specific interference pat- 
tern. The intensity is the modulus- squared of the sum of these three amplitudes plus the 
incoherent addition of /o(1710) and a non-resonant background. The predicted coefficients 
of the /2(1270), 03(1320) and /2(1525) Breit-Wigner amplitudes for an electromagnetic 
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production process are, respectively, +5, -3 and +2 [11,12]. This results in the fit function: 

F(m) = a- |5 £W(/ 2 (1270)) - 3 BW (a° 2 (1320)) + 2 BW(f' 2 (1525))\ 2 + 

+ b-\BW(f (1710))\ 2 + c-U(m), (3) 

where a, b and c are free parameters. 

All the resonance masses and widths were allowed to vary in the fits. The results of the 
fits are shown in Table 1. The quality of both fits, characterized by the x 2 P er number 
of degrees of freedom (see Table 1), is good. However, fit 2 describes the spectrum 
around the i / 2 (1270)/a 2 ) (1320) region better and, unlike fit 1, reproduces the dip between 
/ 2 (1270)/a§(1320) and /2(1525). For this reason and, based on the theoretical motivation, 
fit 2 is preferred and shown in Fig. 2. The background-subtracted mass spectrum is shown 
in Fig. 2b together with the fit. 

The a2(1320) mass in fit 2 is below the PDG value [1]. A similar shift, attributed to 
the destructive interference between / 2 (1270) and ei2(1320), was also seen in a study of 
resonance physics with 77 events [11]. Fit 1 without interference yields a narrow width 
for the combined / 2 (1270)/a2(1320) peak, as also seen by the L3 Collaboration [10]. Fit 
2 with interference yields widths close to the PDG values for all observed resonances. 
The fitted masses for / 2 (1525) and / (1710) are somewhat below the PDG values with 
uncertainties comparable with those of the PDG (Table 1). The quality of a fit without 
the /q(1710) resonance (not shown) yields x 2 / n df= 162/97 and is strongly disfavored. 

The systematic uncertainties of the masses and widths of the resonances, determined 
from the fit shown in Fig. [2], were evaluated by changing the selection cuts and the fitting 
procedure. Variations of minimum track pj>, track pseudorapidity range, track momenta 
by ±0.1%, accepted ti + ix~ mass range around the K s peak and collinearity cuts were done. 
In addition a maximum likelihood fit was used instead of the \ 2 fit an d event selection 
cuts were varied. A check for the possible influence of the J p = + state /o(1500) was 
carried out by including in the fit a Breit-Wigner amplitude of this state interfering with 
the amplitude of the /o(1710). The resulting changes of the fitted values of the mass and 
the width of the / (1710) are included in the systematic uncertainties [13]. The largest 
systematic uncertainties were: fitting with fixed PDG mass and width on /2(1525) affects 
the /o(1710) width by -19 MeV and the largest effect of varying the track momenta on 
the /o(1710) width is +7 MeV. The combined systematic uncertainties are included in 
Table 1. 

The number of events in the Jo (1710) resonance given by the fit is 4058 ± 820, which 
has a 5 standard-deviation statistical significance. This is one of the best /o(1710) signals 
reported. This state is considered to be a glueball candidate [2] . However, if it is the same 
as seen in 77 — > K S K S [9, 10], it is unlikely to be a pure glueball state, since photons 
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can couple in partonic level only to charged objects. Figure [3] compares the results of 
this analysis with other measurements from collider and fixed-target experiments. The 
/o(1710) mass as deduced from the quarkonium decays by the BES Collaboration is signif- 
icantly higher than the values given by all other experiments, including older J/ip-dec&y 
analyses [1]. 

6 Conclusions 

In conclusion, Kg Kg final states were studied in ep collisions at HERA with the ZEUS 
detector. Three enhancements which correspond to /2(1270)/a2(1320), / 2 (1525) and 
/o(1710) were observed. No state heavier than the /o(1710) was observed. The states 
were fitted taking into account the interference pattern predicted by SU(3) symmetry ar- 
guments. The measured masses of the / 2 (1525) and / (1710) states are somewhat below 
the world average, however, the widths are consistent with the PDG values. The /o(1710) 
state, which has a mass consistent with a J PC = ++ glueball candidate, is observed with 
a 5 standard-deviation statistical significance. However, if this state is the same as that 
seen in 77 — > K S K S , it is unlikely to be a pure glueball state. 
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Fit 


No interference 


Interference 


PDG 2007 Values 


xViidf 


96/95 


86/97 


in MeV 


Mass 


Width 


Mass 


Width 


Mass 


Width 


/ 2 (1270) 


1304 ± 6 


61 ± 11 


1268 ± 10 


176 ± 17 


1275.4 ± 1.1 


185.2t|i 


a!j(1320) 


1257 ±9 


114 ± 14 


1318.3 ±0.6 


107 ±5 


72(1525) 


1523 ± 3jJ 


71 ±5±\ 7 


1512 d= 3±J;| 


83 ± 9tl 


1525 ± 5 


73 + _l 


/o(1710) 


1692 ± 6±§ 


125 ± 12t\l 


1701 ± 5±l 


100 ± 24+^ 2 


1724 ± 7 


137 ±8 



Table 1: The measured masses and widths for the / 2 (1270), ag(1320), / 2 (1525) 
and /q(1710) states using K^Kg decays as determined by one fit neglecting interfer- 
ence and another one with interference as predicted by SU(3) symmetry arguments 
included. Both statistical and systematic uncertainties are quoted. The systematic 
uncertainty for the / 2 (1270) /a° (1320) peak is expected to be significant and it is not 
listed. Also quoted are the PDG values for comparison. 
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Figure 1: The measured iv + tt~ invariant-mass distribution for events with at least 
two Kg candidates (dots). The shaded area represents the signal window used for 
K^Kg pair reconstruction. The fit performed (see text) is displayed as a solid line. 
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Figure 2: (a) The measured KgK® invariant-mass spectrum (dots). The solid 
line is the result of the fit described as fit 2 in the text (Eq.(3)) and the dashed line 
represents the background function, (b) Background- subtracted KgKg invariant- 
mass spectrum (dots); the result of the fit is shown as a solid line. 
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Figure 3: Comparison of the present mass and width measurements of the 
/ 2 (1525) and /o(1710) states with other selected measurements [lj. The bands show 
the PDG values and error estimates. 
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